Equation of state of non-additive rf-dimensional hard-sphere mixtures 
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An equation of state for a multicomponent mixture of non-additive hard spheres in d dimensions 
is proposed. It yields a rather simple density dependence and constitutes a natural extension of 
the equation of state for additive hard spheres proposed by us [A. Santos, S. B. Yuste, and M. 
Lopez de Haro, Mol. Phys. 96, 1 (1999)]. The proposal relies on the known exact second and third 
virial coefficients and requires as input the compressibility factor of the one-component system. A 
comparison is carried out both to another recent theoretical proposal based on a similar philosophy 
and to the available exact results and simulation data in d = 1, 2, and 3. Good general agreement 
with the reported values of the virial coefficients and of the compressibility factor of binary mixtures 
is observed, especially for high asymmetries and/or positive nonadditivities. 



I. INTRODUCTION AND A BRIEF REVIEW 
OF THE LITERATURE 

The structure of a dense fluid is known to be largely 
determined by the repulsive intcrmolccular forces, so it is 
not surprising that hard-core potentials have been exten- 
sively employed to model simple fluids and fluid mixtures. 
A noteviforthy aspect of these models is the fact that in 
some instances both exact and approximate analytical 
results may be derived for the structural and thermody- 
namic properties, which in turn serve as a starting point 
for the treatment of more sophisticated or complex mod- 
els. 

Certainly a vast majority of the published work on 
hard-core (rods, disks, spheres, and hyperspheres) fluid 
mixtures pertains to binary systems and to the so-called 
additive hard-core interaction, namely the one in which 
the distance of closest approach (denoted by aij ) between 
the centers of two interacting particles, one of species 
i and the other of species j, is the arithmetic mean of 
the diameters of both particles Ci and CTj , respectively. 
Apart from the initial impetus that took place in the 
60s, recently interest in this kind of systems (in particular 
mixtures of hard spheres) has experienced an increasing 
growth in connection with entropy driven phase transi- 
tions and the demixing problem. On the other hand, 
non-additive hard-core mixtures, where the distance of 
closest approach between particles of different species is 
no longer the arithmetic mean referred to above, have 
received less attention, in spite of their in principle more 
versatility to deal with interesting aspects occurring in 
real systems (such as liquid-vapor equilibrium or fluid- 
fluid phase separation) and of their potential use as ref- 
erence systems in perturbation calculations on the ther- 
modynamic and structural properties of, say, Lennard- 
Jones mixtures. Nevertheless, the study of non-additive 
systems goes back fifty yearaiS*^ and is still a rapidly 
developing and challenging problem. 

As mentioned in the paper by Ballone et al.,"^ 
where the relevant references may be found, experimen- 
tal work on alloys, aqueous electrolyte solutions, and 
molten salts suggests that hetero-coordination and homo- 
coordination may be interpreted in terms of excluded vol- 



ume effects due to non-additivity of the repulsive part 
of the intermolecular potential. In particular, positive 
non-additivity leads naturally to demixing in hard-sphere 
mixtures, so that some of the experimental findings of 
phase separation in the above mentioned (real) systems 
may be accounted for by using a model of a binary mix- 
ture of (positive) non-additive hard spheres. On the 
other hand, negative non-additivity seems to account well 
for chemical short-range order in amorphous and liquid 
binary mixtures with preferred hetero-coordination.'^ 

On the theoretical side, the first exact result on the 
equation of state (EOS) for a non-additive mixture is 
that of a binary mixture of hard rods (d = 1) restricted 
to nearest- neighbor interactions. Although it is usually 
atributed to Lebowitz and Zomick,^ it was already im- 
plicit in earlier work by Prigogine and Lafleur^ and by 
Kikuchif^ and even Lebowitz and Zomick point out that 
the thermodynamic functions of this system appear in the 
thesis presented in 1966 by C. C. Carter (cf. Ref. 9 in@). 
Very recently, Corti and Bowles have rederived this result 
in an appendix of a paper,'' where they also provide ex- 
act geometrical relationships for non-additive mixtures 
(see also an alternative rederivation in Prof. Penrose's 
webpage§.). It is also worth mentioning that in the paper 
by Kikuchij^ a proof is given that no phase transition 
may occur in a onc-dimcnsional binary mixture irrespec- 
tive of the form of the interaction potential, provided it is 
unbounded. The opposite limit of high spatial dimension 
has been considered by Carmesin et al.,^ who showed that 
at sufficiently high density and with positive nonadditiv- 
ity, a binary mixture of non-additive hard hyperspheres 
decomposes into two coexisting phases. 

A very popular model of a non-additive binary mixture 
with positive non-additivity was introduced by Widom 
and Rowlinson in 1970."'^° This model is equivalent to a 
one-component penetrable sphere model. In the sym- 
metric version of the model, referred to as the Widom- 
Rowlinson (WR) model, one has ai — a2 ~ Q and 
CTi2 > 0. With this simplification Widom and Rowl- 
inson derived exactly the EOS in the one-dimensional 
case, where it predicts no phase transition. For d — Z 
the model was solved in the mean-field approximation. 
In the same paper, but for the case of high asymmetry 
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(i.e., when ai ^ 0, (T2 = 0, and d = 3), Widom and 
Rowlinson also determined an approximate condition for 
the spinodal curve. It is interesting to point out that this 
case of high asymmetry corresponds with the Asakura- 
Oosawa model, ^ often used to discuss polymer colloid 
mixtures and where the notion of a depletion potential 
was introduced. This model and refinements of the same 
have received a lot of attention (including fairly recent 
work) in connection with the demixing problem and the 
question of effective potentials.— 

The impact of the WR model cannot be overempha- 
sized as it has motivated a great amount of later work. 
A rigorous proof that a phase transition may exist in the 
WR model in d = 2 was provided by Ruelle,^^ who also 
indicated that a similar procedure may be followed to 
prove the existence of a phase transition in the WR model 
in ci = 3 and higher dimensions. Frisch and Carlier— per- 
formed molecular dynamics simulation for a hard-square 
mixture in the WR limit and concluded that it presented 
a first order phase transition. Melnyck et ali^ obtained 
the first ten virial coefficients of the WR hard-sphere mix- 
ture in the Percus-Yevick (PY) approximation (the first 
five of which are exact), while Straley et alm^ computed 
the virial coefficients of the WR model for oriented hard 
squares and hard cubes. Widom and Stillinger''^& gener- 
alized the scaled particle theory (SPT) for a pure fluid-i 
to the case of the WR model in an arbitrary dimensional- 
ity and Guerrero et ali^ exploited the equivalence of the 
penetrable sphere model and the WR model to obtain 
the direct and total correlation functions for the model 
where the Mayer function is a Gaussian and for the hard- 
sphere interaction in the mean field, PY, and hypernetted 
chain approximations. In the cases of d — 1 and d = 3, 
the WR model was solved in the PY approximation by 
Ahn and LebowitZflS. while the SPT was considered by 
Bergmanns22i The latter theory for the WR model in two 
dimensions was addressed in an appendix of the paper 
by Tenne and Bergmann^Si in which they examined the 
SPT for a non-additive hard-disk binary mixture. Trans- 
port properties for the WR hard-sphere binary mixture 
were computed by Karkheck and Stelli22i Later, Borgelt 
et ali^ and Luo et al.^ performed simulations on the 
hard-sphere WR mixture and found better agreement 
with mean-field results than with PY results. More re- 
cently, the same model has been the subject of investi- 
gations related to its universality class to the location 
of the critical point and the computation of the coexis- 
tence curve r& to the development of an integral equation 
theory that includes the first few terms in the density ex- 
pansion of the direct correlation function into the closure 
approximation,^^ to the (partial) total and direct corre- 
lation functions^ through accurate Monte Carlo simula- 
tions, and to the static and dynamic behavior near the 
consolute critical point obtained from molecular dynam- 
ics simulations.^^ 

A theoretical approach that has been extensively used 
in connection with non-additive hard-core mixtures is the 
SPT. Apart from the papers quoted above, Bergmann^^ 



has considered the SPT in one dimension and compared 
it with the exact solution. Tenne and Bergmann exam- 
ined the SPT for d = 3 both for positive non-additivitj*^ 
(where they computed the critical density and the critical 
non-additivity) and negative non-additivity.'^^ Bearman 
and Mazo also considered the SPT for a symmetric bi- 
nary mixture of non-additive hard-disks'^^ and pointed 
out that the phase transition predicted by Tenne and 
Bergmann in Ref.|^for negative non-additivity was spu- 
rious. The same authors^ introduced a simpler version 
of the SPT for d = 2 and d = 3 which is consistent with 
the SPT of additive mixtures in the appropriate limit 
but still presents some other difficulties. Some of these 
difficulties were addressed by Schaink)2^ who introduced 
an EOS for a binary mixture valid for small values of the 
non-additivity. A comparison of SPT predictions and 
simulation data may be found in Ehrenberg et alm^ 

The use of computer simulation, both molecular dy- 
namics (MD) and Monte Carlo (MC), as well as of the 
usual integral equation approach of liquid state theory 
or the perturbation theory (taking either a one compo- 
nent system or a binary additive hard-core mixture as 
the reference system) have also contributed to the inves- 
tigation of the properties of non-additive hard-core mix- 
tures. In the same paper where they presented the ex- 
act solution for the one-dimensional mixture, Lebowitz 
and Zomicls^ also gave the exact solution to the PY 
equation in d = 1 and a partial solution to the PY 
equation in the three-dimensional case. A mathemat- 
ical analysis of these two solutions was later given by 
Penrose and Lebowitz. Perry and Silbert-- also gave 
an approximate solution to the PY equation in d = 
3 which confirmed the earlier results of Lebowitz and 
Zomick. For equimolar and symmetric hard-sphere mix- 
tures with negative non-additivity, Nixon and Silbert^ 
solved the PY equation, which they found to improve 
its agreement with simulation data as the negative non- 
additivity increased. Equimolar symmetric binary mix- 
tures have been studied by Gazzillo^-""''^^'''^ He has con- 
sidered the PY approximation^ and also other closures 
(the Martynov-Sarkisov,--' the Ballone-Pastore-Galli- 
Gazzillc'* and the modified Verleli^ii^ closures). In Ref. 

he also addressed a ternary mixture with negative 
non-additivity that had been studied earlier through MD 
simulation by Schaink,^'*» while he and his collaborators^ 
were apparently the first to obtain simulation (MC) data 
for an asymmetric hard-sphere binary system. In study- 
ing binary non-additive Lennard- Jones mixtures using 
the RHNC approximation, Anta and Kahl'*^ obtained 
the non-additive hard-sphere bridge functions by solv- 
ing the corresponding PY equation. Lomba et alm^ 
used a generalized modified Verlet closure to study fluid- 
fluid phase separation in symmetric non-additive hard- 
sphere mixtures, obtaining good agreement with their 
own MC simulation data for the phase diagram. Kahl 
et alm^ studied a variety of symmetric binary mixtures 
of non-additive hard spheres (both with positive and 
negative non-additivity) by solving the Ornstein-Zernike 
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equation with a modified hypernetted-chain-type clo- 
sure. Recently, Sierra and DudaiS, considered the PY and 
Martynov-Sarkisov closures to study symmetric mixtures 
of non-additive hard spheres adsorbed on a disordered 
hard-sphere matrix, while Duda et al.*^ used MC simula- 
tions to study fluid-fluid phase equilibria and interfacial 
properties of non-additive binary hard-sphere mixtures 
adsorbed in a slit pore. The structure and the thermody- 
namics of non-additive hard-sphere mixtures under con- 
finement have also been the subject of a recent study by 
PcUicane et al.^ who used both integral equations and 
computer simulations. 

Melnyck and Sawford^- reported MD simulation data 
on a symmetric binary non-additive hard-sphere mixture 
with positive non-additivity and using perturbation the- 
ory derived an EOS for this kind of systems which they 
named MIXl. Such EOS was later extended to cope 
with asymmetric mixtures by Schaink and Hoheisel^ 
At about the same time as the Melnyck and Sawford cal- 
culations, Adams and McDonald^^ performed MC simu- 
lations on binary symmetric hard-sphere mixtures with 
negative non-additivity. Later on, Dickinson^ performed 
MD simulations on two equimolar non-additive binary 
hard-disk mixtures. In 1989, Amar— computed the coex- 
istence curve for the system studied in Ref . |^ using MC 
simulation. Hoheisel^ studied a symmetric equimolar 
binary mixture of non-additive soft spheres with (high) 
positive non-additivity through MD and determined the 
critical density. Mountain and Harvey^ conducted both 
MD and MC simulations on binary mixtures of non- 
additive soft disks to study fiuid-fluid coexistence. Ro- 
vere and PastoreSS, extended the work of Ref. [5^ and 
obtained the coexistence curve of an asymmetric binary 
non-additive hard-sphere mixture through MC simula- 
tion. Extensive MC computations on symmetric non- 
additive hard-sphere binary mixtures have been provided 
by Jung et aLr^i^^^ who have derived from them rea- 
sonably accurate (semi-empirical) equations of state for 
these systems. Density functional theory has also been 
applied^ to the computation of the excess free energy of 
an equimolar mixture of non-additive hard disks. Finally, 
recently Hamad has reported MD calculations for asym- 
metric non-additive binary hard-sphere mixtures^ and, 
together with some coworkers, also for binary hard-disk 
mixtures^ Fluid-fluid phase separation in a symmetric 
mixture of non-additive hard spheres with positive non- 
additivity and the phase behavior of non-additive hard- 
core mixtures in two dimensions have been recently the 
subject of MC simulations by Saija et alw^ and by Saija 
and Giaquinta,^^ respectively, while Gozdz^^ performed 
MC simulations to derive accurate results for the criti- 
cal packing fraction at a few values of the non-additivity 
parameter in the case of hard spheres. Castaheda-Priego 
et alm^ studied depletion interactions in mixtures of non- 
additive hard disks, Schmidt^ generalized the fundamen- 
tal measure density functional theory of hard spheres to 
binary mixtures of arbitrary positive and moderate nega- 
tive non-additivity, and Fantoni and PastoreSSi performed 



accurate MC simulations to check the local dependency 
assumption of the bridge functions of an equimolar non- 
additive binary hard-sphere mixture. Fairly recently, 
Buhoti^ used a cluster algorithm to simulate and study 
phase separation in symmetric binary mixtures of non- 
additive hard disks and hard spheres for various (large) 
non-additivities including the limiting case of the WR 
model. 

An alternative route to the derivation of the EOS of 
non-additive hard-sphere mixtures that does not require 
neither the SPT, perturbation theory, the solution of 
integral equations or simulation results, relies on the 
knowledge of virial coefficients and on the use of ex- 
act statistical mechanical relationships. The so-called y- 
expansion for hard particle fluids introduced by Barboy 
and Gelbarti^ is a prominent example of this approach. 
In the case of non-additive hard-sphere mixtures, the 
Barboy-Gelbart EOS involves up to the exact third virial 
coefficients, whose analytical expressions are knowni^^ 
On a different path, HamadZ^ has provided a theory for 
obtaining mixture properties from pure species equations 
of state. In the case of non-additive hard-sphere mix- 
tures, he invokes exact results pertaining to the contact 
values of the radial distribution functionS)2^*2Si2i as well 
as the knowledge of the exact second and third virial 
coefficients. He has also presented a similar approach 
for hard-disk mixtures in Refs. Is^tM A noteworthy as- 
pect of Hamad's proposal is that, due to his use of the 
one-component radial distribution function as a starting 
point, it is geared essentially towards mixtures not very 
asymmetric in size. This proposal has been very recently 
used in connection with the development of a perturba- 
tion theory for fused sphere hard-chain fluids. 

RecentlyjSMi we have proposed an EOS for a mul- 
ticomponent mixture of additive hard-core particles in 
d dimensions. This proposal shares with Hamad's 
approachSli24i2& two aspects. On the one hand, it is ex- 
pressed in terms of the pure species EOS. And on the 
other it starts with a sensible ansatz on the functional 
form of the contact values of the radial distribution func- 
tions. The aim of this paper is to complement Hamad's 
approach in two different veins. The first one concerns 
dimensionality. Here we want to derive an EOS for a 
non-additive hard-core mixture of an arbitrary number 
of components and for any value of d. The second one 
has to do with the fact that when the non-additivity 
parameter vanishes we also want to recover our former 
proposal^*^ for additive multicomponent hard-core mix- 
tures. Our main concern is to try to keep a reasonable 
compromise between the simplicity of the proposal and 
its ability to deal also with highly asymmetric mixtures. 

The paper is organized as follows. In Sec. ^1 we pro- 
vide general expressions for a multicomponent mixture of 
non-additive hard-spheres in d dimensions and some key 
background material (third virial coefficients, for which 
a simple expression for arbitrary dimensionality is pro- 
posed) for the later development. The exact solution in 
the case of a one-dimensional binary mixture as well as 
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other interesting features of this system are presented 
in Appendix Section 11111 contains a brief account of 
Hamad's proposalS^'iS'^^ for the contact values of the ra- 
dial distribution functions and for the compressibility fac- 
tor of the mixture. This is followed in Sec. lIVI bv our own 
proposal, which shares with Hamad's a few features: the 
construction of the EOS via the contact values of the ra- 
dial distribution functions, the dependence of the latter 
on the EOS of the one-component fluid, and the use of 
the third virial coefficients. The results pertaining to spe- 
cial limiting cases are given in Appendix ^ The analysis 
of the fourth, flfth, and sixth virial coefficients and of the 
compressibility factors in one, two, and three dimensions 
is carried out in Sec. The paper is closed in Sec. IVll 
with further discussion and some concluding remarks. 



yields the virial expansion of Z . namely 

oo 

Z(p) = l + ^(i;rfp)"5„+i 

n=l 

N N 

= 1 + VdP ^ BijXiXj + [vdpf ^ Bijk 

N 

+ {vdpf ^ BijMXiXjXkXt + 0{p^). (3) 

Note that, for further convenience, we have introduced 
the coefficients i?„ = w^''"~^^i?„ where i?„ are the usual 
virial coefficients. The composition-independent second, 
third, and fourth (barred) virial coefficients are given by 



II. THIRD VIRIAL COEFFICIENTS 

A. General equations. 

Let us consider an A^-component mixture of hard 
spheres in d dimensions. The hard core of the interaction 
between a sphere of species i and a sphere of species j is 
CTjj . The diameter of a sphere of species i is a a — Ui. In 
general, fiy — ^{ai + crj)(l -I- Ay ), where Ay > — 1 is a 
symmetric matrix with zero diagonal elements (A^^ = 0) 
that characterizes the degree of non-additivity of the in- 
teractions. In the case of a binary mixture (A^ = 2), 
the only non-additivity parameter is A = A12 — A21. 
The compressibility factor of the mixture Z = p/pksT, 
where p is the total number density, p is the pressure, T is 
the temperature, and is the Boltzmann constant, can 
be exactly expressed in terms of the radial distribution 
functions at contact gij as 

N 

Z{p,{xk},{(Jke}) = l + 2'^~\dpJ2 ^'^J'^i 

^9i3{P,{xk},Wke}), (1) 



2-^-14, (4) 

Bijk = — ^ {ck-ijcrfj + Cj-ikcrfk + Ci-jk<Jjk) , (5) 

Bijke — g (ckeiijC^ij + Cj£-ikcr^k + Cie-jk<^jk 

+Cjk,iicrfi + Cik,:ji><y% + Cij^kicri/,) ■ (6) 

This connection between the virial coefficients of the 
mixture and the c's of the density expansion of the con- 
tact values of the radial distribution functions may be 
profitably used to devise sensible approximations. 

For subsequent use in Sees. IIIII and IIVI it is conve- 
nient to consider the special case of a one-component 
fluid {(Tij = a) of packing fraction y — Vdpc'^ ■ In such a 
case, Eqs. ^ and © become 

OQ 

= l + ^6„+iy", (7) 

71=1 

where bn — i?„/cr'-"^^^'' are the (reduced) virial coeffi- 
cients of the one-component hard-sphere fluid. In partic- 
ular, 62 = 2'*-!. 



where Xi — pi/p is the mole fraction of species i, pi is 
the partial number density of particles of species i, and 
Vd = (7r/4)'^/^/r(H-(i/2) is the volume of a d-dimensional 
sphere of unit diameter. Although no general expression 
is known for gij{p, {xk}, {cTki}) = 9i]{p), it can be ex- 
panded in a power series in density as 

N N 

9ij{p) = 

XkCkiij +0{p'). 

k=l k,t=l 

(2) 

The coefficients Ck;ij, Cki-ij, ...are independent of the 
composition of the mixture, but they are in general com- 
plicated nonlinear functions of the diameters cry , aik, 
(^jk, cTki, ■ ■ .Insertion of the expansion |2Jl into Eq. 



B. The one-dimensional case. 

It is worth recalling that, as mentioned in the Intro- 
duction, in the case of a binary {N = 2) one-dimensional 
(d — 1) mixture with nearest-neighbor interactions only 
[which implies that 2cri2 > max(cri, cr2)]) the exact com- 
pressibility factor is knowniiiiSiLS In Appendix ^ we 
provide a summary of the exact solution as well as some 
interesting properties of the same that, to our knowledge, 
have not been reported before. In particular, the coeffi- 
cients Cfe;y for d — I are 

Cl;ll=Cri, C2;ll = 2(712 - (Tl, Cl;l2=0-1. (8) 

The remaining coefficients are obtained from ^ by the 
exchange 1^2. 
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(a) 



(b) 




FIG. 1: (a) Three spheres of species i,j, and k in an 
aligned configuration. The smallest possible separation be- 
tween spheres i and j is + (Jjk- (b) When sphere k 
is removed, the smallest distance between i and j is Oij. 
Thus crfe;ij = (Tifc + Gjk — Oij represents an eff'ective diam- 
eter of sphere k as seen from the point of view of the pair 
ij. In the sketch we have assumed for simplicity that the 
non-additivities are positive. 



The three-dimensional case. 



In three dimensions, the first two terms of the exact 
density expansion of gij are knowni2& After a few simple 
manipulations one may derive from them the result 



Cfe-i 



3 ^"k-Aj 



(9) 



where 



7]k 



(10) 



and it is understood that (Jk-^j > for all sets ijk. 

cr, , and, in case of additive 



D. A generalization to d-dimensions. 



It is tempting to extend Eqs. 
sions as 



and © to d dimen- 



Cfc;, 



"k-ij 



b2 



'k:i] 



(11) 



More specifically, for a binary mixture Eq. yields 



C2;ll - (2ai2 - (Ti)'^ +1^-1 'Ti(2f7i2 - (Ti)'*^l 



Cl;12 = 



1 (2f7i2 -ai)(Tf/CTi2. 



(12) 

Obviously, Eq. lO reduces to Eq. © for d = 1 (62 
63 = 1), while Eq. ifTT)) reduces to Eq. JHl for d = 3 
(&2 = 4, 63 = 10). 

All of the above results are restricted to the situ- 
ation <7k;ij > for any choice of i, j, and fc, i.e., 
2(Ti2 > max((Ti, (T2) in the binary case. This excludes the 
possibility of dealing with mixtures with extremely high 
negative non-additivity in which one sphere of species k 
might 'fit in' between two spheres of species i and j in 
contact. Since for d = 3 and N = 2 the coefficients Ck;ij 
are also known for such mixtures^ we may extend our 
proposal to deal with these cases. If = 2, one has 
specifically 



C2;ll 



Clearly, ai-ij = CTi, aj-ij 

hard spheres, (Jk-ij — crt- Note also that the quantities 
o'k;ij may be given a simple geometrical interpretation. 
Assume that we have three spheres of species i, j, and k 
aligned in the sequence ikj. In such a case, the distance 
of closest approach between the centers of spheres i and 
j is aik + o'jk ■ If the sphere of species k were not there, 
that distance would of course be cry. Therefore ak-ij as 
given by Eq. HIOII represents a kind of effective diameter 
of sphere k, as seen from the point of view of the interac- Fis) = 
tion between spheres i and j. A schematic representation 
of this interpretation is provided in Fig. ^ 



Cl;12 = (2ai2 - ^2)'^ + (I " 1) ^2fTfM2, 

where we have defined 

CT2 = niax(2(Ti2 - (7i,0) . (14) 

With such an extension, we recover the exact values of 
Ck-ij for a binary mixture of hard spheres {d — 3), even 
if CTi > 2cti2 or (72 > 2(712. We emphasize that Eqs. 
((nil-l|III) for d ^ 1 and d 7^ 3 are new. 

E. The two-dimensional case 

While Eq. I|13|) is exact for d = 1 and ci = 3, it is 
only approximate for d ^ 2. For that dimensionality, 
the exact result has been derived by Al-Naafa et alm^ 
After some algebra (and the correction of some typos), 
the coefficients Ck;ij can be written as 



ci;ii = -^CT?, C2;ii = T'^i^ (^12/0-1) , 

Cl;12 = ^(TfG{<7l2/cTl) , 



(15) 



where 63 = ^ - ^ ~ 3.1280 and the functions F{s) 
and G{s) are given by 



^ (4s2 cos 



-1 J_ 

2s 



s > 
< S < 



1 

2 ' 
1 

2 ' 



(16) 



4 [2.,s^ - 2 (2,s^ - 1) COS- J, 



G{s) = <^ -V4s2 - IJ , 



J.„2 
b3* ' 



Some special values of F{s) and 0(5) are 
F(1) = G(1) = 1, 



F(l/2) = 0, G(l/2) 



2 



0<5<i. 
(17 

(18 
(19 



03 




lim 



03 



lim G(s) = — . 



(20 



For a symmetric mixture ((Ti ~ (72), the value s = 
ci2/o'i = 1 corresponds to the one-component casf 
s = 012! iJ\ = ^ corresponds to the threshold valu 
of negative non-additivity (i.e., 1a\i = (Ti = (T2 
A = ~\)^ 9,nd the limit s — cr\2/cri ~* 00 represent 
an infinitely large positive non-additivity (WR model). 

Equation (|13|l with d — 2 can be recast into the forr 
()15|l . except that the functions F{s) and G(s) are ap 
proximated by 



F (,)_/i(2--l)(4- + &3-4), 



< S < i, 



Gapp(s) 



2s ) 



S > 

< s < 



(21 



(22 



This approximation verifies the properties (|19|I - H2Q|I . ex 
cept that now linis^oo Gapp(s) = 2(63 — l)/63, which i 
about 6% higher than the exact value. Figure |2] show 
that Eqs. H21|l and H22|l constitute an excellent approx: 
mation to the exact expressions (|16f) and p7|l . especiall 
for small or moderate values of s. 

The third virial coefficients Bijk for a two-dimension£ 
binary mixture as given by Eq. |SJ| may be cast into th 
form 



Bill = 2ci;ii(7i = bsaf, 



-B112 = - (C2;ll0-i + 2ci;i2Cri2 

= bacrfS (0-12/0-1) , 



(23 



(24 



and similar expressions for ^222 and B221, obtained fror 
the former by the exchange of indices 1 and 2. Here, 

Bis) ^ ^Fis) + ^s'Gis) 



^ [4ns^ - 8s^ - 1) cos- ± 
- (2s2 + 1) V4s2 - IJ , 



16 „4 
363 ' 



S > 

< s < 



FIG. 2: Plot of the functions F{s), G(s), and B(s). The 
solid lines are the exact functions 11611 . 117II . and I25II . while 
the dashed lines are our approximations l|21|l . 122^ . and 126^ . 

Note that B{s) and i?app(s) are practically indistinguishable. 



Using Eqs. H21() and H22() our approximation yields for 
_Bapp(s) the polynomial 



(25) 



-Sapp(s) = 



363 



< s < i. 

(26) 

As also seen in Fig. |21 Eq. H2t)|) is practically indistin- 
guishable from the exact B{s), so that the (small) dis- 
crepancies in Fapp(s) and Gapp(s) with respect to the 
actual F{s) and G(s) almost entirely compensate. There- 
fore, it seems that it is not unreasonable to use Eqs. Hll|l - 
(fHI for all d. 



III. HAMAD'S PROPOSAL FOR THE 
EQUATION OF STATE 

Our goal is to derive an (approximate) EOS for a mul- 
ticomponent mixture of d-dimensional non-additive hard 
spheres. Clearly, this may be achieved if values for the 
gij are provided. But before we engage in this task, let 
us recall in this section a previous simple proposal by 
Hamad. 

HamadS^*2Si2i has proposed a simple and accurate ap- 
proximation for the contact values of the radial distribu- 
tion functions which takes the same form in both d = 2 
and d ~ 3. Generalized to arbitrary dimensionality d and 
in the notation of this paper it reads 



(27) 



Vdp{cr'^) is the packing fraction of the mixture. 



EN r 
z=i Xicrl 



Here, ; 
with [a^- 

By construction, the aproximation H27|l is correct to 
first order in the density (third virial coefhcient). Insert- 
ing the approximation (|27|) into Eq. (Q, we obtain the 
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(generalized) d-dimensional Hamad EOS 



1 



V 



pure 



inXrj) - 1 



(28) 



So far, the one-component function Zpui.c(2/) remains 
free. It should be emphasized that, except for d — 2 
and d — 3, the EOS given by Eq. (|28|1 has been neither 
introduced nor used before. 

The Helmholtz free energy per particle of a mixture, 
a(p), is given by 



l+Y,x,ln{p,Xf)+ ^ [^(p,) _ 1] , (29) 



where \i is the thermal de Broglie wavelength of species 
i. According to Hamad's approximation (|28|) . 

(30) 

where ap^j.(,(2/) is the excess Helmholtz free energy per 
particle of the pure fluid. 



IV. OUR PROPOSAL 

In 1999 we proposed an EOS for a multicomponent 
mixture of additive hard spheres in d-dimensionsS'^ which 
was based on an ansatz related to the contact values of 
the radial distribution functions. One may express this 
ansatz as 



SYH 



ip) = 



1 



5puro(f?) 



1 



(31) 



where 



(additive spheres) (32) 



is a parameter that is independent of density but depends 
on the composition and diameters of the mixture. 

The idea is now to generalize the ansatz given by Eq. 
(|31|l to the non-additive case. As the simplest possible 
extension, we keep the structure of Eq. (|31|) but deter- 
mine the parameters Zij as to reproduce Eq. |2Jl to first 
order in the density. The result is readily found to be 



-1- 



- 1 



(33) 



The following relationship between and Xij exists 



X,, 



03 



(34) 



The ansatz (|31ll supplemented by Eq. (|33|l is, by con- 
struction, accurate for densities low enough as to justify 
the linear approximation gtj w 1 + VdpJ2k ^kCk-ij ■ 
the other hand, the limitations of this truncated expan- 
sion for moderate and large densities are compensated by 
the use of ^purc- Of course, gij = gpurc in the special case 
where all the diameters are identical (aki = cr) , since then 
Ck;ij = (^3/^2)0''' and Zij = 1. All these comments apply 
to Hamad's prescription H27(l as well. On the other hand, 
Eq. H31|l is consistent, but Eq. H27|l is not, with the case 
of an additive mixture in which one of the species, say 
j = 1, is made of point particles, so that gn — (1 — 77)"^. 

When Eqs. I|3H) and (|33|l are inserted into Eq. |^ one 
gets 



1 



ri (63 - 

M - 1] 



(35) 



Equation H35|l is the main result of this paper. As in 
Eq. (|28|l . the EOS of the mixture is expressed in terms 
of that of the one-component system. On the other 
hand, the density dependence in the EOS (|35|) is simpler: 
Z{p) — 1 is expressed as a linear combination of 77/(1 — 77) 
and ^purc(?7) — 1, with coefficients such that the second 
and third virial coefhcients are reproduced. Again, Eq. 
()35|l is accurate for sufficiently low densities, while the 
limitations of the truncated expansion for moderate and 
large densities are compensated by the use of the EOS of 
the pure fluid. 

In the approximation (|35|l . the Helmholtz free energy 
per particle is 



-,SYH 



(P) 



keT 



-l + ^xan(p.Af) -ln(l-r7) 

i 

53((T'*)^2-52S3 



pure 



)B2 



(63-62)(a'')2 



keT ib3-b2){a^) 



d\2 ■ 



(36) 



In principle, to compute B3, one should use the ex- 
act coefficients Ck-ij- However, since to the best of our 
knowledge they are only known for d < 3 and we want 
our proposal to be explicit for any d, we can make use of 
our approximation for them, Eq. I|ll(l . Therefore, with 
this proviso we get 



63 
02 



Efc Xkcrtii<yi-jkcrj-tk 



- 1 



(37) 



In the additive case {ak;ij (Jk), Eq. reduces to Eq. 
(|32|l . Note that both for d = 1 and d — 3 there is no 
difference in the resulting compressibility factor because 
Eq. Hll|) yields the exact result. On the other hand, for 
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other d, use of Eq. (|37|l also leads to Eq. H35|l , but with 
an approximate rather than the exact value for the third 
virial coefficient. 



V. RESULTS 

Once we have derived our approximation for the EOS 
of the mixture, Eq. H35() . it is interesting to examine 
its performance. And since Hamad has carried out a 
comparison between his proposal, Eq. H28|) . and previous 
Qjjgg^64j74j76j77 finding in general that it performs better, 
we will concentrate here on comparing the results ob- 
tained either through Hamad's prescription or through 
ours (in this regard see also Appendix IB|| . Such compar- 
ison seems in order in view of the fact that both proposals 
share many aspects such as the construction of the EOS 
via the contact values of the radial distribution functions, 
its dependence on the EOS of the one-component fluid 
(more specifically on Zp^rc, that remains to be chosen 
freely) and the use of the third virial coefficients. Also, 
and although Hamad's proposal is specific for d — 2 and 
d — 3 and we have extended it to arbitrary d, they main- 
tain the same form in every dimensionality. Specifically, 
we will focus on the fourth and higher virial coefficients 
and on the compressibility factor. To our knowledge, 
and with the exception of the one-dimensional case, in 
which they are known exactly, values of the former are 
rather scarc o^^i?? and refer exclusively to non-additive 
hard spheres (d = 3). 




-0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 
A 

FIG. 3:_Plots of Biii2/cr» = {Tv/6f Bim/a'-' and Bii22/cr» = 
Bii22/u^ versus A for a symmetric three-dimensional 
binary mixture. Circles: exact valuesf— solid lines: Eqs. I42II 
and l|4i-i|l (present approach); dashed hues: Eqs. I|44|l and H45|l 
(Hamad's result). 



A. Fourth and higher virial coefficients 

From Eq. I|35() it is easy to get an approximate expres- 
sion for the nth virial coefficient: 



B 



SYH 



32 ,dyn-3 



bn - h ,d\n-2 



aT-'B2. (38) 



In particular, the composition independent fourth virial 
coefficients are given by 

^-T [afBjki + (JjBikt. + (yt^ijt + (yfBijk) 



rSyh 



4(6: 



64 - h 



{afafBM + afaiB,, + afa^B.^ 



+a^aiBa + afa^B.k + 4^eB,j ) . (39 
In the case of Hamad's approximation, Eq. one has 



In the special case of binary and symmetric [cti = (T2 = 
tj, ai2 = cr(l-l-A)] three-dimensional mixtures, Eqs. H39(l 
and yield 



i?SYH/^« - 64 ( 1 



rSYH 



4A+yA2 
22 



Ia3 

3 



-A ( 10 + 16A + yA^ 



, 48 ^ 22 ^ 2 28 

1 H Ah A^ H 

9 3 9 



--A(l5-l-24A-mA2) 



A-^ 



(42) 



(43) 



i=.w , q , / 9 , 162 , , 144 , . 
(40) Bfn2/^' = 64 ( 1 + -A + — A2 + —A3 



(44) 



bibl 

"66f 



{(^ijCk-ijCi-ij + (jff.Cj.ikCe-ik + crfiCj-iiCk-u 



+'^'jkCi-jkCi-jk + crfgCi-jeCk-je + ^^M^i-M'^j-M) 



(41) 



rH /^9 



= bi\l + 6A 
688 



408 
'2b 
384 



A^ 



672 
'2^ 
256 



A3 



— A" + A^ + A^ , (45) 

25 25 75 / ' ^ ^ 
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0.35 



0.30 -/'' 




FIG. 4: Plot of B4/B2 versus xi for a symmetric three- 
dimensional binary mixture with A = —0.3 and A = 0.3. 
Solid lines: exact values;^'^ dashed line: Eq. I|38^ (present ap- 
proach); dotted line: Eq. 1401 (Hamad's result). 



where 64 — 18.36477 and we have assumed that A > — 
The two coefRcients B1112 and -B1122 have been evalu- 
ated numerically by Saija et al.^"^ Figure compares the 
numerical data for -B1112 and B1122 with the approxima- 
tions H42|l - H45|l . We observe that Hamad's approxima- 
tion for i3iii2 gives an excellent agreement, while ours 
is only qualitatively correct. On the other hand, for 
Si 122 both approximations are inaccurate for large posi- 
tive non-additivities. In any case, BfY-^ is slightly better 
than B}\22 for < A < 0.3. 

Figure ^ shows, also for a symmetric binary mixture 
of non-additive hard spheres, B4/B2 as a function of the 
mole fraction xi for A — —0.3 and A = 0.3, and the cor- 
responding simulation results. We observe that Hamad's 
approximation is better for A = —0.3, while ours is bet- 
ter for A = 0.3. 

As far as we know, the only report of virial coeffi- 
cients beyond the third for the case of an asymmetric 
non-additive hard-sphere mixture is due to Vlasov and 
MastersiS They have computed up to the sixth virial co- 
efficient for a binary mixture of non-additive hard spheres 
of size ratio 0.1 and a positive non-additivity A = 0.1, 
and up to the seventh virial coefficient for a binary (ad- 
ditive) hard-sphere mixture of the same size ratio. In 
Fig. |31 we present a comparison of the results for the 
composition dependence of the ratio of virial coefficients 



1.0 




FIG. 5: Plot of Bn/Bj versus xi (n = 4, 5, 6) for an asym- 
metric three-dimensional binary mixture with aijax =0.1 
and A = (thin lines) and A = 0.1 (thick lines). Solid lines: 
exact values;*^ dashed lines: Eq. 1381 (present approach); dot- 
ted lines: Eq. II40II (Hamad's result). 



Bn/B^^^ [n = 4,5,6) in the case of a binary mix- 
ture of size ratio CT2/(Ti =0.1 and two non-additivities 
(A = 0,0.1) given by Vlasov and Masters^'^ with the 
results that follow from Hamad's prescription and from 
our proposal. The overall superiority of our proposal in 
this case is apparent and more noticeable for the positive 
non-additivity and when n increases. Nevertheless, the 
negative values of the sixth virial coefficient for the small 
region around x = that are obtained with the simula- 
tion, not shown in the figure, are not captured by either 
proposal. 
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FIG. 6; Compressibility factor Z as a function of A for a sym- 
metric mixture of non-additive hard rods with xi = 0.25 at a 
packing fraction -q = 0.5 (upper panel) and for an asymmet- 
ric mixture with xi = 0.25 and a^joY = 2 at 77 = 0.5 (lower 
panel). Solid lines: exact; dashed lines: Eq. 1481 1 (present 
approach); dotted lines: Eq. (1501 (Hamad's result). 



B. Compressibility factor 

Apart from the virial coefficients, the most important 
tests concern the compressibihty factor itself. In view of 
the big number of parameters in these systems, one has to 
make a judicious choice such that the main features of the 



results may be illustrated. In this subsection we provide 
a representative set of data for different dimensionahties 
that will hopefully cater for the above requirement. 



1. Rods (d=l) 

In the case d = 1, one has Zpurc(y) = 1/(1 ~ J/) ^^nd 
bn = 1 for all n, so that our proposal (|35|l is ill-defined. 
To save that singularity and with the aim of preserving 
the scaling property of the exact solution (see Appendix 
let us write 



^purc 



(y) = 



1 



+ e 



y 



h = 1 



(46) 



1 - y VI - 2// 

and set e ^ at the end of the calculations. In that case, 
replacement into Eq. gives 



7SYH 



(P) - 1+ 



V 



1 



1-77 



which, for a binary mixture, becomes 



7SYH 



ip) = 



1 + X1X2 



CTl + 0-2 V 



-A 



(47) 



(48) 



Note that Eq. H48() is equivalent to a series expansion of 
the exact solution in powers of A truncated in the linear 
term. In fact, in view of Eqs. (jA5|) - (|A7|) . it is exact up 
to order O(A^). Also, it is important to point out that 
Eqs. (|47|l and (|48|l hold regardless of the value of e, so 
the limit e — > has not been needed. 

As for Hamad's approximation, we would have 



Xii = 1 



(Ti + (72 



xaA, Xi2 = 1, 



(49) 



and the similar result for X22 obtained from Xn in Eq. 
(|49|l with the usual replacement 1^2. After some alge- 
bra, one finds 



1 + X1X2 



<J\ + (T2 



1+77 



We remark that Eq. (|50|l is exact to first order in A. 



a;icri/((7) 



X2(J2/{<7) 



1 - 77 (1 + 2i±px2A) 1 - ,7 (1 + £i^a:iA) 



A 



(50) 



r 



A comparison of the exact compressibility factor with 
our approximation (|48|l and Hamad's approximation H50j) 
indicates that Eq. (|48|l . being far simpler than Eq. (|5()() . is 
better than the latter for A > 0, both approaches being 
comparably good for A < 0. This is illustrated in Fig.|Hl 
where we display the exact Z as a function of the non- 



additivity parameter for a symmetric (02! (7\ = 1) and 
an asymmetric (CT2/cri — 2) binary mixture of the same 
packing fraction 77 = 0.5, and mole fraction x\ ~ 0.25, 
together with the two theoretical approximations. 
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0.2 0.3 



FIG. 7: Plot of the compressibility factor versus the non- 
additivity parameter A for an equimolar symmetric binary 
mixture of non-additive hard disks at a packing fraction 77 — 
0.4 (upper panel) and for an equimolar asymmetric mixture 
with a^jox = 3 at 77 = 0.4 (lower panel). The solid lines 
are our proposal, Eq. 13511 . and the dashed lines are Hamad's 
proposal, Eq. (I28II . The circles are results from molecular 
dynamics simulationsi— 



2. Disks (d = 2) 

It seems natural to begin with the case of symmetric 
binary mixtures, i.e., mixtures where ui —02-, and to 
investigate the effect of non-additivity. Representative 
results in this respect for an equimolar symmetric binary 
mixture of non-additive hard disks are displayed in the 
upper panncl of Fig. where we have plotted Z as a 
function of the non-additivity parameter A at a packing 
fraction 77 = 0.4. A similar plot of Z versus A is pre- 
sented in the lower panel of Fig. but in this case for an 
equimolar asymmetric mixture = 3) at the same 

packing fraction r\ = 0.4. 

The size ratio dependence of the compressibility factor 
is displayed in Figs. Inland |51 for various combinations of 
mole fraction xi, non-additivity parameter A, and pack- 
ing fraction ry. 

Although in the paper by Al-Naafa et al^^ they eval- 
uated Z^ by taking for Zpure the one that follows from 
our own simple EOS for the hard-disk fluid, in Figs. 




FIG. 8: Plot of the compressibility factor versus the size ratio 
(72/(71 for an equimolar binary mixture of non-additive hard 
disks with A = —0.2 at 77 = 0.6 (upper panel) and for two 
binary additive hard-disk mixtures (A = 0) at 77 = 0.6 (lower 
panel). The solid lines are our proposal, Eq. 1351 . and the 
dashed lines are Hamad's proposal, Eq. I|28|l . The symbols 
are results from molecular dynamics simulations^- 



|3Elwe have considered for both proposals perhaps the 
most accurate EOS available nowadays, namely the one 
due to Luding)2^ 



{1-yY 64{l-yy 



(51) 



Once again we find that the trend observed in c? = 1 
is also present in the case d — 2, namely, that in general 
our proposal performs better than Hamad's, except for 
negative A. It is worth recalling here that Hamad's EOS 
includes the exact third virial coefficient, Eqs. (|23|) - (|25|l . 
while ours makes use of the approximation embodied by 
Eq. 



3. Spheres (d = 3) 

We proceed here as in the case of d = 2. Figure [TUI 
shows Z as a function of A for a symmetric binary mix- 
ture of non-additive hard spheres at the packing fraction 
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FIG. 10: Plot of the compressibility factor versus the non- 
additivity parameter A for a symmetric binary mixture of 
non-additive hard spheres at 77 = 7r/30 and two different com- 
positions. The solid lines are our proposal, Eq. (1351 . and the 
dashed lines are Hamad's proposal, Eq. 12811 . The symbols 
are results from Monte Carlo simulations p—*— 




FIG. 9: Plot of the compressibility factor versus the size ratio 
(T2/cri for three binary mixtures of non-additive hard disks 
with A = 0.2 at t; = 0.4 and x\ — 0.25 (upper panel), xi = 0.5 
(middle panel), and xi = 0.75 (lower panel). The solid lines 
are our proposal, Eq. 13511 . and the dashed lines are Hamad's 
proposal, Eq. 12811 . The circles are results from molecular 
dynamics simulationsi— 



FIG. 11: Plot of the compressibility factor versus the non- 
additivity parameter A for an equimolar asymmetric binary 
mixture of non-additive hard spheres with size ratio (T2/o"i = 
3, at 7; = 0.5. The solid line is our proposal, Eq. 135L and 
the dashed line is Hamad's proposal, Eq. H28|l . The circles are 
results from Monte Carlo simulations,''^ 



77 = 7r/30 ~ 0.105 and for xi — 0.1 and xi — 0.5. Here, as 
in all the rest of the calculations for hard-spheres, ^puic is 
the one corresponding to the Carnahan-Starhng-Kolafa 
(CSK) EOS,^'^ 

CSK,„, . l + ,,,, 

(1 - y) 

In Fig. ^2 we present a plot of Z versus A, but in 



this case for an equimolar asymmetric non-additive hard- 
sphere mixture with a^jox = 3 at the the packing frac- 
tion T] = 0.5. Finally, Fig.^jis a plot of Z as a function 
of the size ratio for different values of xi, A, and density. 
Once more these figures indicate that our proposal in the 
case of d = 3 is superior to Hamad's, save for negative 
non-additivity. 
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FIG. 12: Plot of the compressibility factor versus the size ratio 
aijox for binary mixtures of non-additive hard spheres with 
x\ = 0.5, A — —0.05, T) = 0.5 (upper panel), x\ — 0.25,0.5, 
A = 0.2, 77 = 0.2 (middle panel), and xx = 0.25, 0.5, A = 0.5, 
r] = 0.075 (lower panel) . The solid lines are our proposal, Eq. 
(I35I I. and the dashed lines are Hamad's proposal, Eq. 1281 . 
The symbols are results from Monte Carlo simulations!— 

VI. CONCLUDING REMARKS 

In this paper we have introduced a new proposal for 
the EOS of a multicomponcnt mixture of d-dimensional 
non-additive hard spheres. This proposal is an immediate 
generalization of the one (rather accurate) we developed 
for additive hard spheres to which it immediately reduces 
if the non-additivity parameters are set equal to zero. 
A general prescription for the d-dimensional composition 



independent third virial coefficients of non-additive hard- 
sphere mixtures has also been introduced. It is exact for 
d = 1 and d = 3 and does a very good job also for d = 2. 
In the absence of exact results or simulation data for 
other dimensionalities, its merits in this respect remain 
to be evaluated. 

Our proposal for the EOS involves providing (sensible) 
approximations for the contact values of the radial dis- 
tribution functions that fulfill a few simple requirements. 
On the one hand, they reduce to the pure component 
value 5purc in the appropriate limit and also comply with 
the limit in which one of the species is made of point 
particles that do not occupy volume. On the other hand, 
they yield the exact gij to first order in the density. Op- 
erationally, our proposed EOS for the non-additive mix- 
ture [cf. Eq. H35|l ] is given explicit in terms of the pure 
component EOS, and the second and third virial coeffi- 
cients of the mixture. The former feature is shared with 
other proposals in the literatureiSiZ^iiLS In any case, 
we find that the present EOS does a good job also in the 
non-additive situation, and represents a reasonable com- 
promise between simplicity and accuracy. In comparison 
with Hamad's approach, which is also simple and rea- 
sonably accurate and which we have generalized here to 
arbitrary dimensionality, it has the advantage of being 
able to deal with asymmetric mixtures where the former 
faces greater difficulties. 

Because the full assessment of our proposal involves 
so many facets, there are of course many issues that we 
have not addressed. We have only attempted to illus- 
trate some of the consequences of employing our approx- 
imate EOS. The results in the previous section illustrate 
a trend that we have observed with other values of the 
parameters, namely that in general Hamad's proposal 
does a better job for negative non-addititivities (espe- 
cially as the density is increased) while ours should be 
preferred in the case of positive non-additivities, at least 
for d = 1, d = 2, and d = 3. Nevertheless, one can see 
that the performance of our EOS is reasonably good in 
highly asymmetric mixtures, even for negative A. So in 
some sense, rather than strictly competing, our approach 
and Hamad's are complementary. It is also worth not- 
ing that here we have chosen to take our original recipe 
of the additive casa^S for simplicity, but we could have 
as well considered the more refined ones that we intro- 
duced laterjSi at the expense of more complicated final 
expressions. Also, the choice of 2'pure is free and the re- 
sults of course depend on that choice. Nevertheless, pro- 
vided .Zpurc is reasonably accurate, the qualitative trends 
should not be altered by different choices and this is actu- 
ally the case. For instance, in the analysis of non-additive 
hard disks we took for Zp^ic the one corresponding to 
Luding's EOS.*^ With minor numerical differences, very 
similar results are obtained if Henderson's equationSl or 
our EOSf24 which are both accurate, are used instead. 
Analogously, in the three-dimensional case the results are 
practically the same if the Carnahan-Starling EOS^ is 
used instead of Eq. H52|l . 
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We are fully aware that interesting features such as the 
demixing transition in the case of positive nonadditivity 
(both for symmetric and asymmetric mixtures) remain to 
be dealt with. We expect to examine some of these in the 
future. In any event, irrespective of the illustrative cal- 
culations that we have presented in this paper, we have 
attempted to include a rather comprehensive account of 
previous work on the subject which will hopefully serve 
to provide some perspective and be useful to other re- 
searchers. 
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APPENDIX A: EXACT SOLUTION IN THE 
ONE-DIMENSIONAL BINARY MIXTURE CASE 

In the one-dimensional case {d — 1) with nearest- 
neighbor interactions [which implies that 2cri2 > 
max ((Ti, (T2)], the thermodynamic and structural prop- 
erties of the binary mixture are exactly knowni^&LS 
The EOS relating the density p to the pressure p (in 
units of ksT) and to the diameters cti, 02, and a 12 = 



(fJi -I- cr2)(l -f A) is given by 



1 1 + 42:1x2 (e2«P - 1) - 1 , , 
- = - + " 2^ A + 



(Al) 



where a = (T12 — (cti + cr2)/2 = (cti + (T2)A/2 > 
— min ((Ti/2, (T2/2). Note that if p — > 00 then 77 ^ 1 

for a > 0, while rj ^ [1 — 2|^min(xi, a::2)]~^ for a < 0. 
Equation (jAljl can alternatively be written as 



P = '5 1 2^1: ■ 



1 — 1] \ 1 — rj J 
where ^{xi,w) is the solution to 



(A2) 



<,-^^l + J-^I±ppI^^lE^ (A3) 
or, equivalently. 

In principle, the compressibility factor is a function of 
four parameters: the number density p, the mole fraction 
xi , the size ratio a2 1 o\ , and the non-additivity parameter 



a. However, the scaling relation HA2|) shows that there 
are only two independent parameters: the mole fraction 
x\ and the scaled parameter w = pajiX — r\). More 
specifically [cf. Eq. (|X2|) ]. 



Z{P) 



1—7] \ 1 — T] 



(A5) 



Thus, ^{xi, w) measures the compressibility factor of the 
non-additive mixture, relative to that of an additive mix- 
ture with the same packing fraction. 

The expansion of the scaling function ^(xi, w) in pow- 
ers of w is 



00 

<^>{xi,w) = ^$„(a;i)w;'^ 

n=0 



(A6) 



where the first few terms are 

$0 = 1, $1 = 2a::ia;2, 
$2 = 0, $3 = -4:xlxl, 



f3 = 

8 2^2 



*5 = -^xlxlil - 8X1X2). 



(A7) 



In the limit of very small non-additivity, we can make 
the linear approximation ^(xi,w) fa 1 + 2xiX2W. This 
is a good approximation in the range —0.4 < w < 0.4. 
The asymptotic behaviors of ^{xi,w) are easily derived 
from Eqs. ljA3|l and ljA4p . In the limit w ^ 00 (with 
a > 0), we simply have $(xi, w) 1, while in the limit 
w —00, the result is $(xi,w) —K{xi)w^^ , where 
K{xi) is the solution to 



4xiX2K =1+ ^^1- 4x1X2(1 -e-^i^). (A8) 

Note also that $(2:1, w) is a non-monotonic function of w 
which presents a maximum for a certain value wo{xi) > 
0. 

From Eqs. (|A5|) - (|A7|) it follows that the (exact) second 
and third virial coefficients can be written as 

B2 = {(t) + 2x1X20, Bz = (ct) ((ct) + 4xiX2a) . (A9) 

Further, the fugacity zi = e'^^ (where pi is the chemical 
potential of species 1, again in units of /csT) is given by 
the following expression 



zi = Aipe^if 1 



1 

2xiw$ 



(AlO) 



and a similar expression for Z2- 



1. Absence of phase separation 

Given the values of ui, 0-2, and (T12 (or a), the thermo- 
dynamic state of the mixture is characterized by the pair 
(PI1P2) or, equivalently, by {xi,w = pa/{l — 77)). Here 
we will adopt the latter viewpoint. If there would exist 
phase separation into two distinct phases A and B, the 
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pressure and the chemical potentials should be equal in 
both phases. The pressure condition is equivalent to 



The conditions on the chemical potentials yield 



(All) 



APPENDIX B: SOME SPECIAL LIMITS 

It is interesting to examine the performance of 
Hamad's approximation, Eqs. H27|l and (|28|l . and of our 
proposal, Eqs. (|31|l and (|35|l . in the following special lim- 
its. 



(A12) 



(712 = 



1-a;^ 



1-xf 



These two equations imply 

_ B 
•^1 — ■''1 1 



(A13) 



(A14) 



(A15) 



Given the non-monotonic behavior of as a function of 
w, Eq. I|A15|) has solutions with ^ . However, the 



combination of ljAll|l and (|A15I) means that 



(A16) 



and so the only solution is the trivial one. 



2. Distribution functions at contact 

From Lebowitz and Zomick's paper^ (and after some 
algebra), one can get the contact values 



922 



1 2xiu;$-($-l) 
1 — T] 2x\w ' 



1 2a;2W$-(*-l) 



1-r, 



2x'\w 



312 = \/5ii522e 
Using the expansion HA6|I . one has 
1 



511 = Tl + 2x2W + 2x2^^ + Oiw^)] 



(A17) 



(A18) 



(A19) 



(A20) 



In the limit of extreme negative non-additivity {ai2 = 
or A = —1), one has di = = 0, so that Eq. l(T^ 
yields 

Cl;12 = Ci:22 = C2;ll = C2;12 = 0. 

The above expressions are exact in that limit. Hamad's 
proposal H27|l becomes then 



xiaf 



X2(T2 

7^' 



(B2) 



9n (P) = 5purc im ) , 9^2 {p) = 5purc {m) , 5^2 (P) = 1 > 

where rji = v^px^af is the partial packing fraction of 
species i. Equation ljB3p is the exact result, reflecting 
the fact that in the limit ai2 = the mixture is actually 
made of two mutually independent one-component fluids. 
On the other hand, in our proposal we have 



hxiaf - b2{a'-) 
(63 - h2){(J'^) 



(B4) 



girip) 



1 hzX2ai , 532^10^^-62(0^ 

I_^(53_6,)(^rf)+3P-W (53_52)(^d) ' 

(B5) 



9ir{p) 



I h , . b2 

9puAm-, 



1 - 7^ 63 - 62 



63 - 62 



(B6) 



plus the equivalent expressions obtained by the exchange 
1^2. Equations ljB5|) and (jB6p are only exact to first 
order in the density. 



Widom— Rowlinson limit 



922 = \l + 2xiw + 2xlw^ + 0{w^)] , (A21) 

1 — 7/ 



312 = [1 - 2x1X2^^ -t- 0{i 

1 — 7? '■ 



The WR limit (ui — <72 ^ 0) represents an extreme 
case of positive non-additivity (A — > 00). The coefficients 
Ok\ij are 



(A22) Ci:il = Ci;i2 = C2;12 = C2:22 = 0, Ci;22 = C2;ll = (2(712)''. 

(B7) 
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In this WR limit the packing fraction vanishes, so that 
the relevant density parameter is 77' — v^paf2- In 
Hamad's approximation, 



From (IB 141) , it follows that 



{^'')Xii = ^X2af2, {a^)X,2 = 0, (B8) 



5ii(p) = 5pu 



26; 



-a;2?/ 



51" (P) 



(B9) 



Our approximation yields 

&i 2^2 erf 2 



2(63-62)' 



K)Z12 = 0, 



(BIO) 



^22 



63 - 62 



1 + ^1 -1 
CTl 



(B16) 



(Bll) 



Both approximations differ in gii{p) and .922 (p), but 
these contact values do not contribute to Z{p) in the 
WR limit. The result is in the two cases 



Z{p) = 1 + 2'^xiX2r]', 



(B12) 



which is the mean field result. 

We note that in the one-dimensional case the exact 
result that follows when setting ai = (72 = 0, 77' = pa in 
Eq. l(X5)l is 



Z{p)^<^>{xuv') (d^l). 



Asakura— Oosawa limit 



(B13) 



The Asakura-Oosawa limit consists of setting (T2 ~ 
and (712 = (Ji/2 + R, where R represents the radius of 
gyration. In that case, CTi = (Ti + 2R and CT2 = 2R, so 
that Eq. lO gives 



Z12 



AR/ai 



1 



Further 

2d ' 



2R/(Ji ' 



xiaf, 



(B17) 



xiaf 



B2 = 
and -83 may be 



this limit {a 

-2(1 -xi) ((7i/2 + i?) 

computed from the c^-^ij given in ljB14|i . Therefore, upon 
substitution into Eq. (|35|l . one would get the EOS for the 
Asakura-Oosawa limit. Since the resulting expression is 
not very illuminating, it will be omitted. Similarly, with 
the substitution of Eqs. ||5T5)| - ||BT7|| into Eq. the 
contact values of the radial distribution functions (which 
will be also omitted) follow. The corresponding results 
for this limit in Hamad's proposal are readily derived 
from Eqs. (|34|l and (jB15|l - ljB17|) . and subsequent substi- 
tution into Eqs. and (|2Hl- 

In d = 1, taking the Asakura-Oosawa limit [02 — 0, 
a = i?) in Eq. ljA5|) we have the exact result 



C2;ll = {2RY + 
Cl;22 = (ffl + 2RY, C2;22 = 0, 



Cl;12 = CTl 



cri+2ij' 



lUi(2i?)" 



C2;12 = 0. 



(B14) 



Z{p)^ 



(d=l). 



(B18) 



* "andres@unex.es' 'http : / /www . u nex . es/f isteor/andres/ 1 

^ malopez@servidor.unam.mx 
http: //miquiztli . cie .unam.mx/xml/tc/f t/mlh/ On 
sabbatical leave from Centre de Investigacion en Energi'a, 
U.N.A.M., Temixco, Morelos 62580 (Mexico) ^ 

* santos@unex.es http : //www. unex . es/f isteor/santos/sby 
^ I. Prigogine and S. Lafleur, Bull. Classe Sci. Acad. Roy. 

Belg. 40, 484, 497 (1954). 
^ S. Asakura and F. Oosawa, J. Chem. Phys. 22, 1255 
(1954); J. Polym. Sci. 33, 183 (1958). 



3 R. Kikuchi, J. Chem. Phys. 23, 2327 (1955). 

* P. Ballone, G. Pastore, G. Galli, and D. Gazzillo, Mol. 
Phys. 59, 275 (1986). 

^ D. Gazzillo, G. Pastore, and S. Enzo, J. Phys.: Condens. 

Matter 1, 3469 (1989); D. Gazzillo, G. Pastore, and R. 

Frattini, ibid. 2, 3469 (1990). 
® J. Lebowitz and D. Zomick, J. Chem. Phys. 54, 3335 

(1971). 

D. S. Corti and R. K. Bowles, Mol. Phys. 96, 1623 (1999). 

* O. Penrose, |http://www.ma.hw.ac.uk/'oliver/hardrods/notes.tex| 



17 



We want to thank Prof. Penrose for providing this new 
derivation to us and for a fruitful exchange of e-mails. 
® H.-O. Carmesin, H. L. Frisch, and J. K. Percus, J. Stat. 

Phys. 63, 791 (1991). 
^° B. Widom and J. S. Rowlinson, J. Chem. Phys. 52, 1670 

(1970) . 

" A. Vrij, Pure Appl. Chem. 48, 471 (1976); A. P. Cast, C. 
K. Hall, and W. B. Russel, J. Colloid Interface Sci. 96, 
251 (1983); E. Canessa, M. J. Crimson, and M. Silbert, 
Mol. Phys. 67, 1153 (1989); E. J. Meijer and D. Frenkel, 
Phys. Rev. Lett. 67, 1110 (1991); H. N. W. Lekerkerker, 
W. C.-K. Poon, P. N. Pusey, A. Stroobants, and P. B. War- 
ren, Europhys. Lett. 20, 559 (1992); E. J. Meijer and D. 
Frenkel, J. Chem. Phys. 100, 6873 (1994); A. Vrij, Phys- 
ica A 235, 120 (1997); T. Biben and J. P. Hansen, Physica 
A 235, 142 (1997); M. Dijkstra, Phys. Rev. E 58, 7523 
(1998); M. Dijkstra, J. M. Brader, and R. Evans, J. Phys.: 
Condens. Matter 11, 10079 (1999); A. A. Louis, R. Finken, 
and J.-P. Hansen, Europhys. Lett. 46, 741 (1999); Phys. 
Rev. E 61, R1028 (2000); R. Roth and R. Evans, Euro- 
phys. Lett. 53, 271 (2001); D. Coulding and J.-P. Hansen, 
Mol. Phys. 99, 865 (2001); A. A. Louis and R. Roth, J. 
Phys.: Condens. Matter 13, L777 (2001); R. Roth, R. 
Evans and A. A. Louis, Phys. Rev. E 64, 051202 (2001); 
M. Schmidt and A. R. Denton, ibid. 65, 061410 (2002); T. 
Bickel, J. Chem. Phys. 118, 8960 (2003); R. L. C. Vink 
and J. Horbach, J. Phys.: Condens. Matter 16, S3807 
(2004); R. L. C. Vink, J. Horbach, and K. Binder, arXiv: 
cond-mat/0409099 (2004). 

^2 D. Ruelle, Phys. Rev. Lett. 28, 1019 (1972). 

H. L. Frisch and C. Carlier, Phys. Rev. Lett. 27, 1040 

(1971) . 

T. W. Melnyck, J. L. Rowlinson, and B. L. Sawford, Mol. 
Phys. 24, 809 (1972). 

J. P. Straley, M. A. Cotter, T.-J Lie, and B. Widom, J. 
Chem. Phys. 57, 4484 (1972). 

B. Widom and F. H. StiUinger, J. Chem. Phys. 58, 616 
(1973). 

^"^ H. Reiss, H. L. Frisch, and J. L. Lebowitz, J. Chem. 

Phys. 31, 369 (1959); E. Helfand, H. L. Frisch, and J. L. 

Lebowitz, ibid. 34, 1037 (1961); J. L. Lebowitz, E. Helfand, 

and E. Praestgaard, ihid. 43, 774 (1965). 
^® M. I. Guerrero, J. L. Rowlinson, and B. L. Sawford, Mol. 

Phys. 28, 1603 (1974). 

S. Ahn and J. L. Lebowitz, J. Chem. Phys. 60, 523 (1974). 
^° E. Bergmann, Mol. Phys. 32, 237 (1976). 

R. Tenne and E. Bergmann, J. Chem. Phys. 70, 1952 
(1979). 

J. Karkheck and G. SteU, J. Chem. Phys. 71, 3620 (1979). 
P. Borgelt, C. Hoheisel, and G. SteU, J. Chem. Phys. 92, 
6161 (1990). 

H. Luo, C. Hoheisel, and J. Karkheck, Phys. Rev. A 42, 
4609 (1990). 

G. Johnson, H. Gould, J. Machta, and L. K. Chayes, Phys. 
Rev. Lett. 79, 2612 (1997). 
'^^ C.-Y. Shew and A. Yethiraj, J. Chem. Phys. 104, 7665 
(1996); K. Jagannathan and A. Yethiraj, ibid. 118, 7907 
(2003). 

A. Yethiraj and G. SteU, J. Stat. Phys. 100, 39 (2000). 
R. Fantoni and G. Pastore, Physica A 332, 349 (2004). 
K. Jagannathan and A. Yethiraj, Phys. Rev. Lett. 93, 
015701 (2004). 
^" E. Bergmann, Chem. Phys. Lett. 38, 385 (1976). 

R. Tenne and E. Bergmann, Phys. Rev. A 17, 2036 (1978). 



E. Bergmann and R. Tenne, Chem. Phys. Lett. 56, 310 
(1978). 

R. J. Bearman and R. M. Mazo, J. Chem. Phys. 88, 1235 
(1988); 91, 1227 (1989). 

R. M. Mazo and R. J. Bearman, J. Chem. Phys. 93, 6694 

(1990) . 

H. M. Schaink, Z. Naturf. A 48, 899 (1993). 

V. Ehrenberg, H. M. Schaink, and C. Hoheisel, Physica A 

169, 365 (1990). 

O. Penrose and J. L. Lebowitz, J. Math. Phys. 13, 604 
(1972). 

B. N. Perry and M. Silbert, Mol. Phys. 37, 1823 (1979). 
J. H. Nixon and M. Silbert, Mol. Phys. 52, 207 (1984). 
D. Gazzillo, J. Chem. Phys. 87, 1757 (1987); Mol. Phys. 
64, 535 (1988); D. Gazzillo and G. Pastore, Chem. Phys. 
Lett. 159, 388 (1989). 

D. Gazzillo, J. Chem. Phys. 95, 4565 (1991). 

D. Gazzillo, Mol. Phys. 84, 303 (1995). 

G. A. Martynov and G. N. Sarkisov, Mol. Phys. 49, 1495 
(1983). 

H. M. Schaink, Physica A 210, 113 (1994). 

J. A. Anta and G. Kahl, Mol. Phys. 84, 1273 (1995). 

E. Lomba, M. Alvarez, L. L. Lee, and N. G. Almarza, J. 
Chem. Phys. 104, 4180 (1996). 

G. Kahl, B. Bildstein, and Y. Rosenfeld, Phys. Rev. E 54, 
5391 (1996). 

O. Sierra and Y. Duda, Phys. Lett. A 280, 146 (2001). 
Y. Duda, E. Vakarin, and J. Alejandre, J. Colloid Interface 
Sci. 258, 10 (2003). 

G. Pellicane, C. Caccamo, D. S. Wilson, and L. L. Lee, 
Phys. Rev. E 69, 061202 (2004). 

T. W. Melnyck and B. L. Sawford, Mol. Phys. 29, 891 
(1975). 

H. M. Schaink and C. Hoheisel, J. Chem. Phys. 97, 8561 
(1992). 

D. J. Adams and I. R. McDonald, J. Chem. Phys. 63, 1900 
(1975). 

E. Dickinson, Mol. Phys. 33, 1463 (1977). 
J. G. Amar, Mol. Phys. 67, 739 (1989). 

C. Hoheisel, Phys. Rev. A 41, 2076 (1990). 

R. D. Mountain and A.H. Harvey, J. Chem. Phys. 94, 2238 

(1991) . 

M. Rovere and G. Pastore, J. Phys.: Condens. Matter 6, 
A163 (1994). 

J. Jung, M. S. Jhon, and F. H. Ree, J. Chem. Phys. 100, 
528 (1994). 

J. Jung, M. S. Jhon, and F. H. Ree, J. Chem. Phys. 100, 
9064 (1994). 

J. Jung, M. S. Jhon, and F. H. Ree, J. Chem. Phys. 102, 
1349 (1995). 

Y. Rosenfeld, M. Schmidt, H. Lowen, and P. Tarazona, 
Phys. Rev. E 55, 4245 (1997). 

E. Z. Hamad, Mol. Phys. 91, 371 (1997). 

M. Al-Naafa, J. B. El-Yakubu, and E. Z. Hamad, Fluid 
Phase Equilibria 154, 33 (1999). 

F. Saija, G. Pastore, and P. V. Giaquinta, J. Phys. Chem. 
B 102, 10368 (1998). 

F. Saija and P. V. Giaquinta, J. Chem. Phys. 117, 5780 

(2002) . 

W. T. Gozdz, J. Chem. Phys. 119, 3309 (2003). 

R. Castaneda-Priego, A. Rodn'guez-Lopez, and J. M. 

Mendez-Alcaraz, J. Phys.: Condens. Matter 15, S3393 

(2003) . 

M. Schmidt, J. Phys.: Condens. Matter 16, L351 (2004). 



18 



R. Fantoni and G. Pastore, arXiv: |cond-mat /040249T1 
(2004). 

A. Buhot, arXiv: |cond-mat/04094"T8l (2004). 

B. Barboy and W. M. Gelbart, J. Chem. Phys. 71, 3053 
(1979); J. Stat. Phys. 22, 709 (1980). 

T. Kihara, Rev. Mod. Phys. 27, 412 (1955). 

E. Z. Hamad, J. Chem. Phys. 105, 3229 (1996). 

E. Z. Hamad, J. Chem. Phys. 101, 10195 (1994). 

E. Z. Hamad, J. Chem. Phys. 105, 3222 (1996). 

H. Hammawa and E. Z. Hamad, J. Chem. Soc. Faraday 

Trans. 92, 4943 (1996). 

E. Z. Hamad and G. O. Yahaya, Fluid Phase Equilibria 
168, 59 (2000). 

B. F. Abu-Sharkh, A. Sunaidi, and E. Z. Hamad, J. Chem. 
Phys. 120, 5795 (2004). The idea of using a non-additive 
hard-sphere mixture as the reference system for other mix- 
tures within the perturbation theory of fluids is not new. It 
was already used in the work of Harvey and Prausnitz [A. 
H. Harvey and J. M. Prausnitz, Fluid Phase Equilibria 48, 
1197 (1989)] and in a paper by Kahl [G. Kahl, J. Chem. 
Phys. 93, 5105 (1990)]. 



A. Santos, S. B. Yuste, and M. Lopez de Haro, Mol. Phys. 
96, 1 (1999). 

A. Santos, S. B. Yuste, and M. Lopez de Haro, J. Chem. 
Phys. 117, 5785 (2002). 

F. Saija, G. Fiumara, and P. V. Giaquinta, J. Chem. Phys. 
108, 9098 (1998). 

A. Yu. Vlasov and A. J. Masters, Fluid Phase Equilibria 
212, 183 (2003). 

A. Santos, M. Lopez de Haro, and S. B. Yuste, J. Chem. 
Phys. 103, 4622 (1995); M. Lopez de Haro, A. Santos, and 
S. B. Yuste, Eur. J. Phys. 19, 281 (1998). 
S. Ending, Phys. Rev. E 63, 042201 (2001); S. Ending and 
A. Santos, J. Chem. Phys. 121, 8458 (2004). 
J. Kolafa, unpublished results (1986). This equation of 
state first appeared as equation (4.46) in the review pa- 
per by T. Boublfk and I. Nezbeda, Coll. Czech. Chem. 
Commun. 51, 2301 (1986). 
D. Henderson, Mol. Phys. 30, 971 (1975). 
N. F. Carnahan and K. E. Starling, J. Chem. Phys. 51, 
635 (1969). 



